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Abstract 

Following Frenod and Sonnendriicker ([H]), we consider the finite Larmor radius 
regime for a plasma submitted to a large magnetic field and take into account both 
the quasineutrality and the local thermodynamic equilibrium of the electrons. We 
then rigorously establish the asymptotic gyrokinetic limit of the rescaled and modified 
Vlasov-Poisson system in a three-dimensional setting with the help of an averaging 
lemma. 

Keywords: Gyrokinetic approximation - Vlasov-Poisson equation - Finite Larmor 
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1 Introduction and main results 
1.1 Physical motivation 

We are interested in the behaviour of a plasma {id est a gaz made of ions with individual 
charge Ze and mass mj and electrons with individual charge — e and mass nie, with rrii >> 
me) which is submitted to a large external magnetic field. It is "well-known" that such 
a field induces fast small oscillations for the particles and consequently introduces a new 
small time scale which is very restrictive and inconvenient from the numerical point of 
view. The simulation of such plasmas appears to be primordial since the model can be 
applied to tokamak plasmas from magnetic confinement fusion (like for the ITER project). 

1.1.1 Heuristic study 

Let us give some heuristic formal arguments to investigate the behaviour of the plasma: if 
we consider the motion of one particle (of charge q > 0, mass m, position x and velocity 
v) submitted to an external constant field B, the fundamental principle of mechanics gives 
that: 

dx dv q , /-, 

— =v,— = —(vAB) (1.1) 
dt ' dt ' ^ ' 
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Straightforward calculations show first of all that the parallel velocity, denoted by tiy (that 
is to say the component of the velocity in the direction of the magnetic field) is conserved 
and thanks to the conservation of the kinetic energy, so is the norm of the perpendicular 
velocity v± (the component of the velocity in the perpendicular plane). Actually, we can 
see that the particle moves on a helix whose axis is the direction of the magnetic field. The 
rotation period (around the axis) is the inverse of the cyclotron frequency fi: 

n = \sM (1.2) 

m 

and the radius is the so-called Larmor radius: 

- M ,,3, 

In the case where the magnetic field is very strong, ^1 tends to infinity whereas tends 
to zero. More precisely, if we take \B\ ~ i (with e 0) we have: 

~ e 

The approximation which consists in considering r^, = is the classical guiding center 
approximation (|17)). This means that each particle is assimilated to its "guiding center" 
(in other words its "instantaneous rotation center"), which is equivalent to neglect the very 
fast rotation of the particle around the axis. 

If one also applies some external constant electric field E, a similar computation shows 
that there appears: 

1. an acceleration in the direction of i?. If we consider £^ ~ 1, then: 

^ - 1 (1.4) 

\B\ ^ ' 

2. a drift -^yr in the orthogonal plane. We have: 

E^B 

(1-5) 

This drift, usually called the electric drift is problematic as regards to the issue of plasma 
confinement. It is negligible compared to the acceleration in the direction of B, but in the 
time scale for plasma fusion which is expected to be very long, one can not neglect this 
small drift, since it creates a displacement of order et {t represents the time). 

At last, note also that if the fields are not constant, various other drifts may appear, 
whose order in e is higher than those of the electric drift. 

Actually, the fields considered are neither constant, nor external, but self-induced by 
the plasma itself. The effects we would like to describe are due to the non-linear interaction 
between the particles and the electromagnetic field. 

1.1.2 The mathematical model 

In all the sequel, we assume that the magnetic field is external and constant and we 
suppose that the speed of particles is small compared to the speed of light, so that we 
can use the electrostatic approximation which consists in reducing the Maxwell equations 
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to the Poisson equation. Finally, we decide to opt for a kinetic description for the ions: 
in other words, the time and space scales considered here are such that ions are not at a 
thermodynamic equilibrium and their density is governed by a kinetic equation. 

The basic model usually considered for the ions is the following Vlasov- Poisson system: 

f dtf + v.VJ + {E + vA B).VJ = 

-A^V = f fdv 
ft=o = fo 

where f{t,x,v) is the density of ions, with t £R+,x £R'^ or R'^/Z'^,v E M'^ (usually 
d = 2 or 3), meaning that f{t,x,v)dxdv gives the number of ions in the infinitesimal 
volume [x, x + dx] x [v,v + dv] at time t (note that in this model, electrons are for the 
moment neglected). 

1.1.3 The gyrokinetic approximation 

It is important from a numerical point of view to establish the asymptotic equation when 
\B\ tends to infinity. Indeed, we expect the asymptotic equation to be "easier" to handle: 
only one time and space scale, perhaps less variables in the phase space to deal with... 
The derivation of such equations is usually referred to in the mathematic literature as 
"gyrokinetic approximation". 

Rigorous justifications of these derivations with various time and space observation 
scales have only appeared at the end of the nineties. We refer for instance to the works 
of Brenier ([5]), Frenod and Sonnendriicker (|11|-[T2]). Frenod, Raviart and Sonnendriicker 
([To]), Golse and Saint-Raymond ([IS]- [16]), Saint-Raymond f[2T]-[22]). 

The classical "guiding center approximation" corresponds to the following scaling for 
the Vlasov-Poisson system (from now on and until the end of the paper, i? is a constant 
vector, say for instance B = \ez)'- 

{ dth + v.VJ, + {E, + ^).V,A = 

-A,V, = Jfdv ^'-^^ 

, fe,t=0 = fo 



The articles [111 and |15j show that when e — )• 0, this leads to a one-dimensional kinetic 
equation in the direction of B: 

dtf + v\\.V.,f + E\\.V,f = 

E = -v.y 

-A,V = Jfdv 
ft=o = fo 

Notice that the electric drift does not appear; this was expected since we have seen in 
the formal analysis that this drift was of higher order in e than the other effects. This 
shows in particular that this approximation is not sufficient for the numerical simulation 
of tokamaks. In order to make this drift appear, there exists to our knowledge two main 
possibilities: 

1. one consists in restricting to a 2D problem in the plane orthogonal to B (|15)). 

2. the other consists in rescaling the orthogonal scales in order to get both transport 
and electric drift at the same order (11211. 
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This work directly follows the articles |12) and |10| where the authors considered the 
"finite Larmor radius approximation". This means that the spatial observation scale in 
the plane orthogonal to B is chosen smaller than the one in the parallel direction, more 
precisely with the same order as the Larmor radius r^,, so that one can expect the electric 
drift to appear in the asymptotic equation. 

In some sense, having such a scaling allows the electric field to significantly vary across 
a Larmor radius, which is not the case for instance in (|L6p . Moreover, in this situation, the 
positions of the particles are no longer assimilated to the position of their "guiding center" 
and we will have to perform an average over one fast oscillation period (the so-called 
gyroaverage) in order to get a sort of averaged number density. 

1.2 Scaling and existing results 

The system we are going to study is based on the "finite Larmor radius scaling" and takes 
into account the quasineutrality of the plasma. 

1.2.1 The (refined) mathematical model 

We refer to [12] for a complete discussion on the scaling. Let us recall briefly and quite 
crudely how it works. 

Let L|| be the characteristic length in the direction of the magnetic field and L± be the 
characteristic length in the perpendicular plane. We consider that ~ 1 and L± ~ e and 
define the dimensionless variables xj| = and x'j_ = . In the same fashion we also define 
the dimensionless variables t' and v' with characteristic time and velocities with the same 
order as Ly and introduce the new number density /' defined by ff'(t',x',v') = f{t,x,v) 
(and we define likewise the new electric field and potential EE' (t' , x' ,v') = E{t,x,v) and 
VV'{t' ,x' ,v') = V{t,x,v)). We consider the scaling /,£' ~ 1 and V ^ e. At last, we 
introduce the Debye length of the plasma A/), which appears in the Poisson equation. In 
order to take into account the quasineutrality of the plasma, we take from now on Ad ~ ^/e. 

The Poisson equation states in this scaling: 

- eA,,V: - -^A,,^V: = ^ {nl - <) (1.8) 

where n\ = j f^dv' is the density of ions and the density of electrons. The density 
distribution of ions is normalized so that J f^dy'dx' = 1. 

The main difference between Frenod and Sonnendriicker's model and ours lies in the 
following. Instead of considering a fixed background of electrons, and since ^ << 1, we 
make the usual assumption that the (adiabatic) electrons are instantaneously at a local 
thermodynamic equilibrium, so that their density follows a Boltzmann-Maxwell distribu- 
tion: 

<(x,t)=exp(£^) (1.9) 

where kB is the Boltzmann constant, — e the charge and the temperature of the electrons. 
We consider that , % ~ 1. 

We make the assumption that we are not far from a fixed background of electrons, so 
that we can linearize this expression: 

<(x,t) = l + T/' (1.10) 
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We are obviously aware that this assumption is not reahy satisfactory both from a 
mathematical and physical point of view; let us just consider this as a plain mathematical 
model which allows us to give rigorous justifications in this work. The problem of a fixed 
background of electrons, i.e. = 1, brings actually more interesting formal results; this 
point will be discussed in the last section. 

The Poisson equation can now be written: 

Vi - e^A,, - A,,V^ = I fdv' - I hdv'dx' (1.11) 

The dimensionless system ()1.6p becomes (for the sake of simplicity, we forget the 
primes) : 

{ dtfe + ^.V,A + v^^.VJ, + {E, + ^).V,/, = 

V, - e^A^ii V, - A^^V, = f f,dv - J fodvdx ^ ' > 

. fe,t=0 = fe,0 

with the notation Ax^^ = d^^^ and Axj_ = A — Aa;||, 

the problem being posed for (x_L,X|j,i)) G x T x M'^ (with T = M/Z equipped with 
the restriction of the Lebesgue measure to [0, 1[). 

1.2.2 State of the art about the Finite Larmor Radius Approximation 

Using homogenization arguments, Frenod and Sonnendriicker established the convergence 
in some weak sense of sequences of solutions {fe)e>o of similar systems, in two cases, namely 
in some pseudo 2D case (assuming that nothing depends on and ?;||) and in a 3D case 
when the electric field is external. The main tool used to establish the convergence is the 
"2-scale convergence" introduced by Nguetseng |20| and Allaire [2] that we will recall later 
on. 

1. The 3D case : 

Assume that we deal with an external electric field = E ^ (M. x M.^): 

dtfe + t-ii-V,./, + ^.V,A +{E + ^) .V,/, = 
ft=o = /o 

Frenod and Sonnendriicker proved the following theorem: 

Theorem 1.1. For each e, let be the unique solution of the scaled Vlasov equation 
in Lf°(Ll. Lx y). Then the following convergence holds as e tends to 0; 

f,^ f weak-* LriLl,) (1.13) 
where f £ Lf^L^.^) is the unique solution to: 

dtf + t-ii .V,./ + (^^'" nr)E{t, X + n{-T)v)dT^ .Vxf 
^ (^J^^ R{T)E{t, X + n{-T)v)dT^ .V,/ = 

/|t=0 = ^ (^'" fo{x + n{T)v, R{T)v)dT^ 



5 



denoting by R and TZ the linear operators defined by: 

cos r — sin r 
sin T cos T 
1 



Rir). 

2. The pseudo 2D case: 



7^(r) = (-i?(-7r/2) + R{-7r/2 + r)) 



The Vlasov-Poisson system considered in this case is the following 2D system: 



dtfe + ^.V,./e +Ie, + ^] .V,/e = 



fe\t=0 — fo 

-\/V,,-A^V, = p. 



Pe 



fedv 



(1.14) 

(1.15) 
(1.16) 

(1.17) 



li V = {vx,Vy), f-*" is defined by [vy, —Vx). 

We recall that there exist global weak solutions of Vlasov-Poisson systems in the 
sense of Arsenev (|3j)- 

Assuming here that /o > 0, /o G H Lx,v (for some p > 2) and that the initial 
energy is bounded, Frenod and Sonnendriicker proved the following theorem (we 
voluntarily write an unprecise meta-version of the result) 

Theorem 1.2. For each e, let {f^, E^) be a solution in the sense of Arsenev to ^1.14^ - 

mi). 

Then, up to a subsequence, weakly converges to a function f Moreover, there exists 
a function G such that : 

f= G{t,x + n{T)v,R{T)v)dT (1.18) 
Jo 

and G satisfies : 

dtG + ^ T^ir)£it, T,x + 7^(-r)^>)(^r^ .V,G 

+^ (/ ^(^)'^(*' ^' ^ + T^{-r)v)dT^ .V.,G = 
^1^=0 = fo 

£ = -V<^, -A^> = j G{t,x + n{T)v,R{T)v)dv 
denoting by R and TZ the linear operators defined by : 



Rir) 



cos T — sm T 
sin T cos T 



,7^(T) = (i?(-7r/2) - i?(-7r/2 + r)) 



In this case, we have to introduce an additional variable, the "fast-time" variable r 
which comes from the fact that we need to precisely describe the oscillations in order 
to study the limit in non-linear terms. 

Note that the authors actually developped a generic framework that allows them to deal 
with different scalings and to give a precise approximation at any order. We do not wish 
to do so in our study. 
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1.3 A bit of homogenization theory and some useful definitions 

Let us now precisely state the "2-scale" convergence tools used in this paper. 

Definition. Let X be a separable Banach space, X' be its topological dual space and (., .) 
the duality bracket between X' and X. For all a > 0, denote by Cq(M, X) (respectively 
Lq (M; X') ) the space of a-periodic continuous (respectively ) functions on M with values 
in X. Let q G [1; oo[. 

Given a sequence (ue) of functions belonging to the space L'^ (0, and a function 

U°it,9) G L'i'{0,T;Li{R;X')) we say that 

Ue 2-scale converges to 

if for any function ^ G L''(0, T; Cq(M, X)) we have: 

lim^ (^u,{t),^ (^t,^jj dt = ^ {U^{t,T),-^{t,T))dTdt (1.19) 

Theorem 1.3. Given a sequence (li^) hounded in L^ (0,t;X'), there exists for all a > a 
function C/° G L«'(0, T; (M; X')) such that up to a subsequence, 

Ue 2-scale converges to 

The profile [/^ is called the a-periodic two scale limit of and the link between and 
the weak- * limit u of is given by: 

i / U°dT = u (1.20) 
" Jo 

We also introduce some notations: 

Notations. We define for all p G [1; oo] the space L^ y'. — ij^(T'^, (i7^(M'^))). 

In the same fashion, we define the spaces L^,j., L^^^... 

Let -^277 T space of 2Tr -periodic functions of r which are in L%. 

Let be the space of functions f such that for all infinitely differentiable cut-off 

functions (p G , (pf belongs to Lx- We will say that a sequence (f^) is uniformly bounded 
in L^ if for each compact set K , the sequence of the restrictions to K is uniformly 
bounded in Lx with respect to e (but this bound can depend on K). 

We will also use the same notations for Sobolev spaces W^'^ (s gR). 



1.4 Statement of the result 

In this paper we prove that the 2-scale convergence established in the previous 2D case is 
also true in our 3D framework. The difficulty comes from the fact that there is no uniform 
elliptic regularity for the electric field because of the factor in front of Aa;|| in the Poisson 
equation: 

V, - e^Ax^^V, - Ax^V, = J fdv - j fdvdx 

In particular there is no a priori regularity on xy and therefore no strong compactness. 
Nevertheless, we actually prove that due to the particular form of the asymptotic equation, 
the moments of the solution with respect to v\\ are more regular in xy than the solution 
itself. We can then easily pass to the weak limit. 
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The reason why we have opted for this strange Poisson equation instead of the usual 
one win appear at the end of the next section and especially in the last one. Roughly 
speaking it allows us to "kill" plasma waves which appear in the parallel direction due to 
the quasineutrality. 

Notice that this result is in the same spirit as the proof of the weak stability of the 
Vlasov-Maxwell system by DiPerna and Lions (|8]), where the authors have regularity on 
moments, by opposition to the proof of the weak stability of the Vlasov- Poisson system by 
Arsenev ([3]), where the author has compactness on the electric field. Actually our result 
is a kind of a hybrid one, since we get on one hand regularity with respect to x± by elliptic 
regularity and in the other hand regularity with respect to by averaging. 

We assume here that the initial data (/e,o)e>o satisfy the following conditions: 

• /e,o > (positivity) 

• {fe,o)e>o is uniformly bounded with respect to e in H Lx,v (for some p > 3) and 
for each e, / fefidxdv = 1. 

• The initial energy is uniformly bounded with respect to e: 
j f,^o\v\^dvdx + e J KV^ + e J \'^x^V,^o\'^dx + j |V^||K,oPrfx^ <C 

Theorem 1.4. For each e, let {f^,E^) in Lf{Ll.y n L^,^) x Lf{Ll.) he a global weak 
solution in the sense of Arsenev to jil.l^) . Then up to a subsequence we have the following 
convergence as e tends to 0: 



fefi weakly-* converges to /o G L^. ^ (1-21) 
/, 2- scale converges to F G n J) (1.22) 

2-scale converges to £: G Lf°°(L^^^(L^/2(wif ))) (1.23) 

Moreover, there exists a function G G L'^{L\ ^ n L^^v) such that: 

F{t, T, X, v) = G{t, X + niT)v, R{t)v) (1.24) 

and {G,£) is solution to: 

dtG + .V^G + ^ T^{r)£{t, T,x + n{-T)v)dT^ .V^G 

^14=0 = /o 

£ = {-V^V,0), V-A^V = j G{t,x + n{T)v,R{T)v)dv - j fodvdx 
denoting by R and TZ the linear operators defined by: 



R{t) 



cos T — sin T 
sin r cos r 
1 



, 7^(T) = (i2(-vr/2) - R{-tt/2 + r)) 
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As it has been said, for the proof of this theorem, we will first prove a proposition which 
gives the regularity of moments in v^^ of the solution. For this, we use an averaging lemma. 
The beginning of the proof is very similar to the proof in the 2D case, but we will give it 
again for the sake of completeness. 

Remarks. 1. The assumption on the initial energy may, at first sight, look a bit restric- 
tive but in the "usual" Vlasov-Poisson scaling, it only means that the inital electric 
potential and field are bounded in L^. 

2. The constant q = 3 will come quite naturally from Lemma [2. 31 and Proposition 13.11 

3. This theorem implies that for a given non-negative initial data G|j=q = Gq in ^ n 
Lx,v (with p > 3) and satisfying the energy bound, the asymptotic system admits at 
least one global weak solution G G Lf^(L],_^nLS,«). With the additional assumptions 
on the inital data: 

Go G W^^l 
11(1 + |?;|^)Go||l- < CO, 
||(l + |i;|^)I)Go||L- <oo. 

we are actually able to prove the uniqueness of the solution, using the same ideas 
than Degond in [7] (and also used afterwards by Saint-Raymond in a gyrokinetic 
context (|21))). Hence, it means that if the whole sequence {fe,o) satisfies the same 
additional estimates, uniformly with respect to e, and weak * converges to some /o 
then there is also convergence for the whole sequence (fe)- 



2 A priori uniform estimates for the scaled Vlasov-Poisson 
system 

2.1 Conservation of norms and energy for the scaled system 

In this section we give a priori estimates which are very classical for the Vlasov-Poisson 
system (used for example in |11| . |12| . |15|). In order to recall how one can get them, 
we will give some formal computations. If one wants to have rigorous proofs, one should 
deal with smooth and compactly supported functions, namely with a sequence (/")n>o of 
solutions of some regularized Vlasov-Poisson equations then pass to the limit (that is the 
way one can clasically build a global weak solution in the sense of Arsenev ([3])). 

First, as usual for such Vlasov equations, norms are conserved (we work here at a 
fixed e): 

Lemma 2.1. For all I < p < oo, 

yt>o,\\mu.^^ < wmui^ (2.1) 

Moreover, fo>Oif and only i/Vt > 0,f{t) > (referred to as the maximum principle) 

That precisely means that if /o G Lx,v, then / G Lf[L^^y). 
Let us now compute the energy for the scaled system: 

Lemma 2.2. We have the estimate: 
£^{t) = (^j f,\v\'^dvdx + e J V^dx + e j \V ^ ^V^? dx + j IV^,, Kl^dx^ < f.(0) (2.2) 
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In particular if there exists C > independent of e such that <Se(0) < C, then: 

fe\v\'^dvdx < C (2.3) 



Formal proof . We multiply the scaled Vlasov equation by \v\'^ and integrate with respect 
to X and v. 

j dtfe\v\^dvdx + j E,.Vi,f,\v\'^dvdx = ^ (^j f^\v\'^dvdx^ - 2 j E,{x).vf,dvdx = 

We then integrate the Vlasov equation with respect to v. We get the so called conser- 
vation of charge: 

^ (^J fdv^ + V,,|. (^J fv^^dv^ + y fv^dv^ = (2.4) 

Therefore, we have: 

j E^{x).vf^dvdx = - j {\/x^Ve,eVx^^V^).vfedvdx 

= y K (V,^ . {f,v^) + eV^ii . {f,v\\ )) dvdx 
= —e Vedtfedvdx 



Finally, using the Poisson equation, we get: 

= -e( [ VAV.dx + / V, , VAV^ , V,dx + / V ^AV.V ^.V.dx 



ld_ 

''2di 



-e^^i I V^dx+ I |V,,F,|2(ix + e2 / \V,,V,\^dx 



Thus it comes: 



^ (^j fM'^dvdx + e j V^dx + e j \V^^V,\'^dx + j \V^^^V,\'^dx^ =0 (2.5) 



□ 



2.2 Regularity of the electric field 

Let us recall a classical lemma obtained by a standard real interpolation argument: 
Lemma 2.3. Let f{x,v) be a mesurahle positive function on x M^. Then: 

/ r- N 3/2 y r. N 3/4 

{j f{x,v)dv\ dx < C\\f\\%^^^ U \v\^fdxdv\ 



(2.6) 
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Proof. For any i? > 0, we can write the following decomposition: 



f,dv = I fdv+ fdv 

\v\<R J\v\>R 



< Ci?2||/||^3 + -l / \v\'fdv 
Then we can take R such that i?^||/||i3 = f \v\'^fdv so that we get: 



1/6 / /• \ 1/2 

p3 j„, \ ( / l„,|2 . 



fdv<C^J fdvj ^J\v\'fdvj (2.7) 

We then raise the quantities to the power 3/2, integrate with respect to x and use Holder's 
inequality which gives the estimate. 

□ 

By conservation of the norm and the uniform bound on the intial energy, Lemmas 
\2J\ and O entail that: 



G LriLl/') (2.8) 

and the norm is bounded uniformly with respect to e. 

We now use the Poisson equation to compute the regularity of the electric field. Let us 
recall that: 



V, - e^A^..V, - A^^Ve = Pe- Podx 



Lemma 2.4. With the above notations and assumptions: 

is uniformly bounded with respect to e in L'^ {Lx^^^ {Wxf^'^)) 

Proof. Let e > and t > be fixed. For the sake of simplicity we write V instead of 
and E instead of E^. 

For any function f{x^^,x±), define the rescaled function f{z,x±) by 

so that: 

ll/(^,^±) 11.3/2 = ||/(2;||,xx)|| 3/2 (2.9) 
The Poisson equation becomes: 

V - A,V - Ax^V = p-e-i J Podx 
and the scaled electric field is given by: 

E= [-V,V,-Vx^V 

Since Pe{t, ., .) and are uniformly bounded in Lj , standard results of elliptic regu- 
larity on the torus x show that there exists C > independent of e such that: 



11^11^2 3/2 < C 



p-ei I Podx 



,3/2 
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Remark. Notice here that due to the dilatation of order ^ in the parahel direction, being 
periodic in this direction does not make things easier. 

Thanks to (|2J]) we get: 



1^11^2 3/. < C 



p- / podx 



<Co 



with Co independent of e. 
Consequently, we have: 



Finally from ()2.9p we get 



^||^3/2(^M/2) < \\E\\^l,3/2 < Co 



□ 



We can see as expected that the regularity of the electric field with respect to the 
variable is not sufficient to get some strong compactness. 

Remarks. 1. We can write the identity: 

- A^^V, = -A,^{Id - e^A,,! - A,^)-i (^p, - j p,dx^ (2.10) 

so that, thanks to elliptic estimates on the torus T^, E Lx^^^ {Wxf^"^). Consequently, 

Q /'2 \ 3 /2 

dx^^Ve is bounded in Lxj_ {Wx^^ )• This implies that y = —edx^^V^ tends to zero 
in the sense of distributions. 

2. A typical function ip^ such that ip^^ is bounded in L'^ and ^ip^^ is bounded in W~^'^ is 
the oscillating function cos(^x). This indicates that E^^^i oscillates with a frequency 
of order ^ in the parallel direction. 

3. If we work with the usual Poisson equation 

-e^Ax..Ve - Ax^V^ = Pe- Pedx 



we only get homogeneous estimates for and we have not been able to deal with 
such anisotropic estimates in the following of the paper (namely in the estimates 
of Proposition 13. ip . Roughly speaking, if 1/ is a solution of the Poisson equation 
-AV = p with p G L3/2(m3)^ 

we can only say that V G W'^'^/'^{^^) (the homogeneous 
Sobolev space) and not l^^'^/^. 

4. This difficulty seems to be not only a technical one, but appears to be linked to the 
existence of plasma waves (with frequence and magnitude of order -^) in the parallel 
direction which prevents us from passing directly to the limit e — ?• (see |19] and last 
section) . 



3 Proof of Theorem 11.41 

Proof. The first two steps are identical to the one given in |12) . For the sake of completeness 
we recall here the main arguments and refer to |12] for the details. 
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Step 1: Deriving the constraint equation 



First of all, since (/e) is bounded in L^(L;|, ^CiLx^v), Theorem 11.31 shows that for all a > 0: 
/, 2-scale converges to F„ G L°°(0, T; L^(M; „)) 
Let ^{t,T, x,v) be an a-periodic oscillating test function in r and define: 

= ^{t, -,x,v) 

We start by writing the weak formulation of the scaled Vlasov equation against Since 



we get the following equation: 



0, 



+ ,(y^q,Y ) dtdxdv 



fo^{0,0,x,v)dxdv 



Multiply then by e and pass up to a subsequence to the (2-scale) limit. We get the so 
called constraint equation for the a-periodic profile F^. 



drFa + v±.VxEa + vA e^.V„F„ = 0, 
which means that F^ is constant along the characteristics: 



dV 

Ik 
dX 

d7 



(3.1) 

(3.2) 
(3.3) 



A straightforward calculation therefore shows that there exists F^ G L°^{0,T, Lx,v) 
such that: 

F^{t, T, X, v) = F^{t, X + n{T)v, R{t)v) (3.4) 



with: 



Rir) 



cos r — sin r 
sin T cos r 
1 



,7^(r) 



-R{-tt/2) + R{-tt/2 + t)) 



sin T cos T — 1 
i.e. TZ{t) = 1 — COST sinr . 

0_ 

Since R and TZ are 27r-periodic, we will consider the 27r profile: indeed if a and 27r 
were incommensurable, F^ could not depend on r and consequently we would have no 
information on the oscillations. 
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Step 2: Filtering the essential oscillation 

We now look for the equation satisfied by := G; we introduce the filtered function g^: 

g,{t, X, v) = Mt, X + n{-t/e)v, R{-t/e)v) (3.5) 

(meaning that we have removed the oscillations) 
We easily compute the equation satisfied by g^: 

dt9e + ^^11 ■V.^ge + n{t/e)E,{t, x + n{-t/e)v).V^g, (3.6) 
+R(t/e)E^{t, X + n{-t/e)v).V^g^ = 

Remark. Note here that g^: 2-scale converges to G, and since it does not depend on r, it 
also weakly converges to G. 



Step 3: Getting some regularity on moments 

From now on, the goal is to get some compactness for the moments of with respect to 
I'll . The main tool we have in mind is the following averaging lemma proved by Bezard in 
[6], which is a refined version of the fundamental result of DiPerna, Lions and Meyer ([9]): 

Theorem 3.1. Let 1 < p < 2. Let f,g& L^^dt ig) (g) dv) be solutions of the following 
transport equation 

dtf + v.VJ = (/ - Ai,,)-/2(/ - A^r/^g (3.7) 

with m G M+,r G [0,1[. ThenM'^ G C^(M'^), p^{t,x) = J f{t,x,v)^{v)dv £ W^'f (M x M"') 
where 

(3.8) 



(1 + m)p' 
Moreover, 

WP'i'Wwt'l^iRxM.'i) — C {\\f\\LP{dtmx(S)dv) + \\9\\LP{dt(S>dx^dv)) (3-9) 
(C is a positive constant independent of f and g) 

Averaging lemmas are an important feature of transport equations: since the transport 
equation (13. 7p is hyperbolic, one can obviously not expect the solution / to be more regular 
than the right hand side or the inital data. Nevertheless, if one considers the averaged 
quantity /Jiir, one can actually notice a gain of regularity. This phenomenon was first 
observed independently by Golse, Perthame and Sentis ([H]) and Agoshkov ([T]) then was 
formulated in a precise way for the first time by Golse, Lions, Perthame and Sentis (see 
|13|): it is referred to as "velocity averaging". There exists many refined versions of these 
results and numerous interesting applications in kinetic theory, but we shall not dwell on 
that. We simply point out that this tool has been successfully applied to Vlasov equations, 
for instance to prove the existence of global weak solutions to the Vlasov-Maxwell system 
as it has been done by DiPerna and Lions ([H]). 

These results have been proved for functions with values in M. Here, for our purpose, 
we need a new version of L^ averaging lemma for functions with values in some Sobolev 
space W^'P{M^) [k G N*). We prove the following result, which is sufficient in our case 
(probably an analogous of Bezard's optimal result is also true): 

Lemma 3.1. Let 1 < p < +oo and A G M. Let f,g £ L^^^{Wy'^) be solutions of the 
following transport equation 

dtf + v.V,f = {L-A,r/^g (3.10) 
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with m £ M+. Then G C^iM"^), pis,{t,x) = J f{t,x,v)^{v)dv £ W^'^iWy'^) for any s 
such that 



and 



Moreover, 



s<sp = —^-— forp^2 (3.12) 
(1 + m)p' 



(C is a positive constant independent of f and g) 

Sketch of proof. We prove the result in the stationnary case only: 

v.V,f = {I-A,r/^g (3.14) 

By standard arguments (see |13) ) the general case then follows. 

The following estimate is obvious for q = 1 oi q = +oo (and actually we can not expect 
any smoothing effect) : 



For p = 2, we prove the result as in Golse-Lions-Perthame-Sentis |13| . We denote by ^ 
(resp. rj) the Fourier variable associated to x (resp. y). 
The only point is to notice (using Fubini's inequality): 

MmHy^) = jii + \i?y'^ + :F^,r,f^{v)dv'j drjdi 

The proof is then identical and we get for s = 2(-\Xm) ' 

MnsiHy^) < C (||/||l2_^(ha) + Wghi^iHy^)) (3-16) 
Finally the general case 1 < p < +oo is obtained by complex interpolation [1]. 

□ 

Equipped with this tool, we can now prove that moments in v^^ are more regular with 
respect to t and xy than the solution itself. 

Proposition 3.1. For each e > 0, let be a function in Ll.^ri Lx,v (with p > 3) bounded 
uniformly with respect to e and satisfying: 

dt9e + ^^11 -^xge + n{t/e)E,{t, X + n{-t/e)v).Vxg, 
+R{t/e)E^{t, X + n{-t/e)v).V^g^ = 

with the electric field uniformly bounded in Lf{Lx )■ 
Let ^ G P(M). Define 



Ve{t,x,v±) = J ge{t,x,v)^{vii)dvi\ 

Then, 

?7e is uniformly bounded in W^'J^^^i^^{W~^'^]^j^J (3-17) 
for 7 G]1; 2[ defined ^ = | + ^ C'l^d some s G]0; 1[ (depending on 
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Proof. • The first step is to localize the equation. Let K be the cartesian product of 
compact sets: 

K = [0, T] X K^ii X K^^ X i^^ii X K^^ 

We now consider some positive smooth function ^{t, x^^,x±,v^^,v±) which is and 
which satisfies the condition: 

^' = outside K (3.18) 

Noticing that: 

diVa; (n(t/e)E,{t, X + n{-t/e)v)) + div„ {R{t/e)E^{t, x + n{-t/e)v)) = 
The equation satisfied by is the following one: 

dt{ge^) + u^.V,ig,<^) = - V.,.{n{t/e)E,{t,x + n{-t/e)v)g,^) 

' V ' 

(1) 

- V^.{R{t/e)E,{t, X + TZ{-t/e)v)g,^) -dt{^)g, - t;,, .V^(^>)(7. 

^ V ' 

(2) 

+n{t/e)E,{t, X + n{-t/e)v).V^{<^)g, + R{t/e)E,{t, x + n{-t/e)v).V^{^)g, 

The idea is now to consider this equation as a kinetic equation with respect to the vari- 
ables (ijXijjWii) and with values in an abstract Banach space (which will be Wxl^vj_)- 
We then only study the first two terms of the right-hand side (noticing that the other 
terms have more regularity than these ones). 

From now on, for the sake of simplicity and readability, we will write and W^'^ 
norms without always specifying that they are taken on the compact support of 

• Estimate on the first term (1) 

Since E^ does not depend on v, we have: 

E,eLr{Ll/^{L^{Ll/^'))) 
In particular if we restrict to compact supports: 

The second point is that the differential operator applied in (1) involves only deriva- 
tives with respect to the x± variable and not in the parallel direction: this remark is 
fundamental for using the averaging lemma 13.11 (indeed, the case of a full derivative 
in X|| can not be handled). 

Holder's inequality simply implies that: 

\\n{t/€)E,it,x + ni-t/e)v)gMLZ, „, < \\E,it,x + ni-t/e)v)n ||<7.|Lp 

(3.19) 

where ^ = | + |- Hence: 

\\Vx.(n{t/e)ES,x+n{-t/e)v)g,^)\\^-i., <\\E,{t,x+n{-t/e)v)^\\ 3/2 I^Jl? „ 

(3.20) 
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Notice that the change of variables (x, v) {x + TZ{s)v, v) has unit Jacobian for all 
s E M, so that: 

\\E,{t,x + n{-t/e)v)^^3/2 =\\E,{t,x)\\^y2 (3.21) 



So finally we have, after integrating in and thanks to Holder's inequality: 

II V^. {n{t/e)E,{t, X + 'R{-t/e)v)g,^) ||^, , i,^ . 

t,X II ,U|| V ''^'^ ,^'_[_ / 

< C'l|-Ee(i,x)|| 3/2 ^ 3/2 \\g^\\lV (LP ) 

and C is a constant independent of e. 

Remark. The regularity of (1) with respect to v± is not optimal (since it involves 
no derivative in v±_ for g^^). Nevertheless we are interested in the regularity of the 
whole right hand side, and we will see that the term (2) has this regularity in v_\_. 

Estimate on the second term (2) 

By the same method one gets: 

II V„. {R{t/e)E,{t, X + n{-t/e)v)g,^) \\ 

Finally we see that the right hand side is uniformly bounded in: 

r,X||,loc\ D||,loc\ X ± ,v ± ,loc ■> 

Regularity of the moments 

By lemma \3A\ . for all ^ G C^, the moment: 



r]e{t,x,v±) = J ge{t,x,v)'^{v\'^)dv\\ 
is then uniformly bounded in the space W^'J ioc(^2r^'J loc) ^o^' s > with 



^ < 27' 



XX,f_L,loc/ 

1 



□ 



We can now prove that the sequence of moments r]^ is compact in a space of distributions 
which is the dual of some space where the sequence (Ef^) is uniformly bounded. 

Corollary 3.1. There exists 9 g]0, 1[ and rj G Wf^'^^i^^{W~^'^^^j^J such that for all £, > 0, 

up to a subsequence: 

7?, ^ 77 Strongly m i^?,,oc(^Lioc(W^x^x';f '.D) (3-22) 
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Proof. By assumption on the initial data, there exists q > 3 such that /o G L^^^; thanks to 
the a priori L'^ estimate, we get g^, G L'^{L%^v). Define 7 by: 

1 _ 2 1 

7 3 g 

The previous lemma shows that for some s > 0: 

Ve e W^t';';|ioc(^xl^;Jx,loc) uniformly in e 
Since g^^ E L^^^^{L'i^v) and ^ has compact support, we get by Holder's inequality: 

Since ^ > | > | and | < 5, there exists e]0, 1[ such that 

1 i-e e 

3^7 

By interpolation ([!]) we deduce that: 
This implies that: 

Ve G ^^t:foc(^x||(^-"x,l,loc)) uniformly in e 

Ve G ^?,loc«':ioc( V:i,loc)) uniformly in e 

We then use the following refined interpolation result proved by Simon in [23], which 
is, roughly speaking, an anisotropic adaptation of the classical Riesz-Frechet-Kolmogorov 
criterion for compactness in L*': 

Theorem 3.2. Let 1 < p < 00 and s > 0. Let T > and X, B, Y be three Banach Spaces 
such that X C B C Y and with X compactly embedded in B. Let F be a bounded set of 
Lf([0,T],X) n W/'^([0,r],y). Then F is relatively compact in Lf([0,T],S). 

This entails, thanks to Sobolev's embeddings, that the sequence (tj^) is strongly rela- 
tively compact in Lii,M,,Xoc{W^!^%,)), for aU i > 0. 

□ 

From now on, we consider such that 9 + ^, < 1, which is of course possible since 9 < 1. 

Remark. Following the remark in Step 2 and by uniqueness of the limit in the sense of 
distributions, we get: 

V = j G^{vii)dvii 
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Step 4: Passing to the weak limit 

We will first need a technical lemma which is obtained directly from the 2-scale convergence 
of E,. 

Lemma 3.2. Up to a subsequence, 

• TZ{t/e)E^{t,x + TZ{—t/ e)v) 2-scale converges to ^ J^^ Tl{T)£{t,T,x + Tl{—T)v)dT G 

• R{t/e)E^{t,x + lZ{—t/e)v) 2-scale converges to j^^ R{T)£{t,T,x + TZ{—T)v)dT E 
LT{Lf^,AL%^{wlf'^)) 

Proof. E^ is uniformly bounded in L'l" [L^^^^ {Wxf^'^)) , so there exists £ G L1° {Lx^^^ {Wxf^'^) 
such that Ef: 2 scale converge to £. 

We take ^{t,T,x) a 27r-periodic w.r.t. r test function and use the 2 scale convergence 
ofE,: 



n{t/e)E^{t, x + 7^(-^/e)^;).^'(^, t/e, x)dtdx 
E^{t, x)}n{t/e)^{t, t/e, X - n{-t / e)v)dtdx 

^[ [ £{t,T,x)}TZ{T)^{t,T,x -n{-T)v)dtdTdx 







2-K 



— / / £{t,T,x + n{-T)v)}TZ{T)'^{t,T,x)dtdTdx 
2vr J Jq 



The proof is the same for R{t/e)E^{t,x + TZ{—t/e)v). 



□ 



Now, we can write the weak formulation of the kinetic equation (|3.6p against a smooth 
test function of the form ^[t,x,v±_)^{v\\/) with compact support. If we can pass to the 
limit for such test functions, then by density it will be also the case for all test functions. 

Noticing that div^; v\\ = and that 

div^ (n{t/e)E,{t, X + n{-t/e)v)) + div„ {R{t/e)E,{t, x + n{-t/e)v)) = 
we get: 

J X, v^)^{v\^)) + .Vxim + n{t/e)E,{t, X + n{-t/e)v).Vx{m 

+R{t/e)E^{t,x + n{-t/e)v)SJ^{<^^!)^gedtdx^dx\\dv^dv\\ 

uq^{Q, X, v±)"^{v\i)dxdv 



We can easily take weak limits in the linear part dtge + "^11 xQt 



19 



Consider now the "non linear" term: 

TZ{t/e)E^{t,x + n{-t/e)v)SJ.^^{t,x,v±)ge^{v\\)dtdx^dx\\dv±dv\\ = 
TZ{t/e)E^{t,x + TZ{-t/e)v).\7x^{t,x,v±) (^j g^^{v\\)dv\^ dtdx±dx^^dv± 

The convergence of this term can be estabhshed by the strong/weak convergence prin- 
ciple. Nevertheless, we have to carefully use this technique to get the result and we will 
explicitly evaluate the difference: 



n{t/e)E^{t, X + n{-t/e)v).V^<^{t, x, v ^)riedtdx ^dx\\dv ^ 

- f — f n{T)£{t,T,x + n{-T)v)dT.Vx<^'ndtdxdvj_ 

< J {n{t/e)E^{t, X + n{-t/e)v) - ^ J T^{-r)£{t, r, x + 7^(-r)^;)(ir^ .V ^^r]dtdxdv i_ 

+ j TZ{t/e)E^{t,x + TZ{-t/e)v).V^^{t,x,v±){r]^-r])dtdx_i_dxiidv± 

The first term of the right hand side converges to zero because of the 2-scale convergence 
of E^ (Lemma 13. 2p . We can control the second term by: 

(these norms are actually taken on the compact support of ^ but we do not write it for 
the sake of simplicity) 

Using the fact that E^ is uniformly bounded in L^j^^(Lx(j^(wf^^'^^^)) (this is an easy 
consequence of Lemma 12. 4p and Corollary 13.1^ 

we can deduce that 

j 'R{t/e)E^{t,x + TZ{-t/e)v).V^^{t,x,vx) [i]^ - if) dtdx_i_dx^\dvj_ 

The proof is of course the same for the other non-linear term: 

R{t/e)E,{t, X + 7^(-^/e)^;).V„5, 

To conclude let us compute the asymptotic equation satisfied by the 2-scale limit of 
denoted by V. We take ^{t,T,x) a 27r-periodic w.r.t. r test function. We write the weak 
formulation of the Poisson equation: 

V^^^V,V:,^^^{t,t/€,x)dtdx + j V^^V,V^^^{t,t/€,x)dtdx 
fe{t,x,v)"^{t,t/e,x)dtdvdx — / ( / fodvdx] "^{t,t/e,x)dtdvdx 
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We then pass to the 2 scale hmit: 



— f f V{t,T,x)Va:^'^{t,T,x)dTdtdx + 
2vr J Jq 

+ TTZ i i Vx^V{t,T,x)Vx^'^{t,T,x)dTdtdx = — / / F{t,T, x,v)'^{t,T, x)dtdvdx 

2tt 



2vr 7 7o ■ ■ ^ ' ■ ' ' ■ ^ ' • ' ' - 2tt J JO 
1 

2^ 



1 

2^ 



27r 



y J fodvdx^ ^{t,T,x)dTdvdx 

— f f G{t,T,x + TZ{T)v,R{T)v)'^{t,T,x)dTdvdx 
2vr J Jo 

fodydx^ ^{t,T, x)dTdvdx 
from which we get the "Poisson" equation given in Theorem 11.41 
V - A±V = J G{t, X + n{T)v, R{T)v)dv - j h 



dvdx 





Moreover since E^^± = — V^^^V^ and thanks to Remark 1 fohowing Lemma I2.4| we 
easily get if we pass to the two-scale limit: 

□ 

4 Concluding comments 

4.1 Comments on the result 

Finally we can see as in |12) (namely by performing the change of variables x = Xc — v-^ 
and looking at the new equations in the so-called gyro-variables {t,Xc,v)) that the drift 
involving the electric field in the asymptotic "kinetic" equation corresponds to the electric 
drift that we mentioned in the introduction and which was expected to appear. We also 
notice that the Poisson equation we get in the asymptotic system is the same than the one 
used in the numerical simulations of tokamak plasmas (see for example the GYSELA code in 
|18|). Nevertheless, physicists do not get it in the same formal way: they claim that it only 
expresses the quasineutrality of the plasma (there is no "real" Poisson equation involved) 
and the perpendicular laplacian happens to appear due from the so-called "polarization 
drift" (|18|. see also |24) for a physical reference on the subject). It would be interesting to 
justify such a computation from a mathematical point of view. 

At last, we wish to point out a really unpleasant feature of our model, which is that 
there is no parallel dynamics. 

4.2 An alternative model 

Let us give some comments on the gyrokinetic approximation of the system (|4.2p which 
consists in considering a population of electrons in a fixed background of ions: 

ni = [ fodxdv (4.1) 
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Actually, quite surprisingly, this model engenders more interesting physical properties: 
in this case the parallel component of the electric field does not vanish but appears as 
a pressure in the end (which may bring difficulties both in the study of the asymptotic 
system and in numerical simulations). 



-e^A^^l^e - A^^F, = / fdv - J fodvdx 

. fe,t=0 = fe,0 

With the same computations as the present paper, we get: 

/e 2-scale converges to F 
E^ 2-scale converges to £ 

In a formal sense, there exists a function G such that: 

F{t, T, X, v) = G{t, X + n{T)v, R{t)v) 

and {G,£) is solution to: 

dtG + v\\ .V^G + ^ T^{r)8{t, r, x + 7^(-r)^;)(^r^ .V^G 

^14=0 = /o 

£ = {-yi_V,£\\), -/\i.V = j G{t,x + n{T)v,R{T)v)dv- j h 
still denoting by R and IZ the linear operators defined by: 



(4.2) 



(4.3) 
(4.4) 



(4.5) 



dvdx 



R{t) 



cos T — sin T 
sin r cos r 
1 



,7^(T) = (i?(-7r/2)-i?(-7r/2 + T)) 



The parallel component has to be seen as a pressure (or the Lagrange multiplier) asso- 
ciated to the "incompressibility" constraint fj2 f G{t,x + TZ{T)v,R{T)v)dvdx± = J fodvdx 

Let us just give a few words on the difficulties that arise with this model. The Poisson 
equation can be restated as: 



e2A.,„K-A.,,K 



so that thanks to the linearity of the Poisson equation we can study separately two equa- 
tions. The first one states: 



dvdx I 



(4.7) 



For this part of the electric potential we get the same estimates as in lemma [231 Indeed, 
J (J fedv — f fedvdx±) dx± = so that we can use elliptic estimates on the torus T^^. 
Consequently this electric potential does not give birth to any parallel dynamics, like in 
Theorem 11.41 
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The second one is: 



- e^A^iiK^ = j Uvdx^ - J fodvdx (4.8) 

^ V ' 

only depends on 

This equation, associated to the one giving the electric field E^^y = —edx^^V^^, is similar 
to the one studied by Grenier in |19) , coupled to a Vlasov equation describing a quasineutral 
plasma. In this case it was shown that there exist plasma waves with both temporal and 
spatial oscillations with frequency and magnitude of order Because of these waves, 
it is much more difficult to pass to the limit in order to get a kinetic equation. Grenier 
managed to prove the convergence only for distribution functions with special form and got 
in the end a Euler-like system with an electric field interpreted as a Lagrange multiplier. 
Hence, in our case, this part of the electric field may engender an non-zero . 

For these reasons, it seems much harder to expect to prove a rigorous result similar to 
Theorem 11.41 with such a model. 



4.3 Prospects 

A way to pass to the limit in this latest case would be to use a relative entropy method 
like in the papers of Brenier ([5]) and Golse and Saint-Raymond (|16)). This will be the 
object of a forthcoming work. 

Another interesting issue would be to consider a "true" Boltzmann-Maxwell distribution 
for the electrons (not linearized like in this paper) and perform an asymptotic analysis, 
maybe also with a relative entropy method. 

Acknowledements. The author would like to thank L. Saint-Raymond for having given 
him the idea to write this paper and for her careful reading and advice. The author is also 
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